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Abstract 

In e + e~ collisions at centre-of-mass energies around 91 GeV, D*° mesons have been recon- 
structed using data collected with the OPAL detector at LEP. The hadronisation fraction has 
been measured to be 

/(c^D*°) = 0.218 ±0.054 ±0.045 ±0.007, 

where the errors correspond to the statistical and systematic errors specific to this analysis, 
and to systematic uncertainties from externally measured branching fractions, respectively. 
Together with previous OPAL measurements of the hadronisation fractions of other charmed 
mesons, this value is used to investigate the relative production of observed vector and pseu- 
doscalar charmed mesons in Z° — > cc decays. The production ratio is determined to be 

pf = v/(V+P) = 0.57 ± 0.05 . 

The relative primary production of vector and pseudoscalar mesons, Py im , is studied in the 
context of the production and decay of orbitally excited charmed resonances. The first mea- 
surement of the inclusive D* + production rate in hadronic Z° decays is presented. 
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1 Introduction 



A meson with no orbital angular momentum can be a vector (V) state with spin 1 or a pseu- 
doscalar (P) state with spin 0. Their relative production Py=V/(V+P) is sensitive to non- 
perturbative effects in the hadronisation process and cannot be calculated exactly. However, 
several models |IJ have been proposed which predict this ratio. A simple spin counting picture, 
where the abundance of a particular state is proportional to its number of spin degrees of free- 
domQ, predicts a value of Py=0.75 when only vector and pseudoscalar meson production are 
considered. More sophisticated models take into account the mass difference between vector 
and pseudoscalar mesons, as well as the masses of the constituent quarks. In general, they 
predict Py to be less than 0.75. 

The situation is complicated by the presence of mesons with non-zero orbital angular mo- 
mentum. The observed Py ratio thus includes those ground state mesons which have been 
produced in decays of the excited states. Experimentally, this effective ratio P v s is easier to 
measure. If the production rates of the excited states are known, the primary ratio Py im , 
corrected for any effects from excited states, can be calculated. The LEP accelerator, where 
numerous Z° — > qq decays have been observed between 1989 and 1995, provides a facility to 
study the Py ratio. Values around Py ff ~0.75 and P v s £z0.55 have been measured for B meson 
production p| and in the charm sector || [|], respectively. In both cases, a model dependent 
evaluation suggests the values of p^ nm and P{f to be similar ||. For light mesons, P v s has 
been estimated to be between 0.4 and 0.5 0. This paper focuses on a study of both P v s and 
Py im for charmed meson production. 

In the charm system, measurements of D°, D + , and D* + production in Z°— >cc decays have 
been used in the previous determinations of Pf @, §]. The D*° meson has so far not been 
observed in Z° decays, since it only decays via the emission of a photon or a 7r° meson, which 
are difficult to reconstruct experimentally. Therefore, the determinations of Py values have so 
far relied on the assumption of isospin invariance, which suggests equal D*° and D* + production 
rates. 

In this paper, a first measurement of the hadronisation fraction /(c — » D*°) in Z° — > cc 
decays is presented. The analysis is based on more than 4 million hadronic Z° decays recorded 
with the OPAL detector at the LEP accelerator in the years 1990 to 1995. The D*° mesons 
are reconstructed in the decay channels D*° — > D°7 and D*° — > D°7r°. The same techniques 
are applied in the D+7 final state for the reconstruction of D* + mesons, which at LEP have 
only been observed in leptonic D+ decays @. The /(c^D*°) measurement is used to test the 
assumption of isospin invariance. Together with previously published OPAL measurements of 
the other charmed non-strange pseudoscalar and vector mesons || 0, a value of P{f is derived. 



A recent OPAL measurement of the production of excited charmed mesons |10[ is used to 
investigate the ratio P v Tim . A model independent formula for P v nm is derived, and the validity 
of a simple spin counting model for the fragmentation process is tested. 

The paper is organized as follows: Section!^ contains a brief description of the OPAL detector 



1 "Spin counting" in this paper refers to a model where the relative primary production of mesons with the 
same quark content and the same orbital angular momentum is given according to the corresponding numbers 
of spin degrees of freedom. 
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and the event selection. Section |3] describes the reconstruction of D*° decays, the selection of 
signal candidates, and the background determination. In section [|, the measurement of the 
hadronisation fraction /(c— >D*°) and its systematic uncertainties are discussed, while section [5] 
describes the measurement of D* + production. Section |6| contains a determination of both P v s 
and Py im for charmed meson production in Z° — ■> cc decays and an interpretation of the results. 



2 The OPAL Detector and Event Selection 



A complete description of the OPAL detector is given elsewhere [nj. Here, only the components 



of importance for this analysis are reviewed. Tracking of charged particles is performed by a 
silicon microvertex detector, a vertex detector, a jet chamber, and a set of drift chambers that 
measure the coordinate of tracks along the direction of the beam linef] (^-chambers), positioned 
inside a solenoid that provides a uniform magnetic field of 0.435 T parallel to the beam direction. 

The barrel electromagnetic calorimeter, which covers the polar angle range of | cos#| < 0.82, 
is mounted outside the magnet coil. It consists of a cylindrical array of 9,440 lead glass blocks of 
24.6 radiation lengths thickness pointing approximately to the interaction region. The overall 
energy resolution is improved by correcting for the energy lost in showers initiated in the 
material in front of the calorimeter. Such showers are detected by thin gas detectors (presampler 
detectors) situated in front of the lead glass blocks, and by time-of-flight scintillators located 
between the presampler and the magnet coil in the polar angle range | cos#| < 0.72. The regions 
0.82 < | cos ^| < 0.98 are covered by the endcap electromagnetic calorimeters with lead glass 
blocks oriented parallel to the beam direction. The magnet return yoke is instrumented as a 
hadron calorimeter. Four layers of muon chambers are mounted outside the hadron calorimeter. 

The criteria for selecting hadronic Z° decays are based on reconstructed tracks in the central 



detector and on the energy distribution in the calorimeter |12| . Charged particle tracks need 
to have at least 20 jet chamber hits, a momentum component in the xy plane of at least 
0.15 GeV, a total momentum of less than 65 GeV, and a distance of closest approach to the 
beam axis of less than 5 cm. The hadronic Z° event selection efficiency of these requirements 
is (98.7 ± 0.4)% 0. Of the events recorded with the OPAL detector between 1990 and 1995, 
4.32 million satisfy the event selection criteria. The primary vertex is reconstructed from the 
charged tracks in the event and constrained with the known average beam position and the xy 
width of the e + e~ collision point. 

Samples of simulated hadronic events are used for the determination of selection criteria and 
for the calculation of selection efficiencies. They have been generated using the JETSET 7.4 
Monte Carlo model ]13| with parameters tuned to reproduce the OPAL data [14 . The frag- 



mentation of heavy quarks is parametrised by the fragmentation function of Peterson et al. [JT5 



The simulated events are then processed by the detector simulation program [161 and by the 
same reconstruction algorithm which is also applied to the data. 



2 The OPAL coordinate system is denned with positive z along the electron beam direction and the x axis 
horizontal; 9 and </> & r e the polar and azimuthal angles, respectively. 
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3 Reconstruction and Selection of D*° Candidates 



In Z° decays, D*° mesons are dominantly produced in Z° — > cc events and in the decay of bottom 
hadrons. The D*° mesons are reconstructed in the decay modes D*° — > D°7 and D*° — > D°7r°. The 
transition photon or tt° meson is expected to be dominantly produced in the core of the jet that 
contains the D*°. Thus, for the reconstruction of D*°— >D°7 decays, a considerable background 
level is expected from photons originating in decays of 7r° mesons from other sources, while 
the reconstruction of D*° — > D°7r° decays suffers from combinatorial background in the 7r° — > 77 
reconstruction. The reconstruction proceeds by first finding D° candidates, which are then 
combined with a photon or a 7r°. After a loose preselection, a likelihood method is used for the 
final selection of candidates. In this section, the selection and the background determination 
are described. 



3.1 Reconstruction of D° Mesons, Photons, and 7r° Mesons 

Both the preselection for the D°, photon, and ir° reconstruction and the likelihood for the D*° 
selection are based on the variables which are explained below. 

The decay mode D°— >K~7r + is used to reconstruct D° mesons. Initially, all possible pairs of 
oppositely charged tracks are formed, assigning kaon and pion masses to the tracks. Variables 
for the selection of D° candidates are: 

• the invariant mass m D o and the scaled energy x D o = E D o / E\, cam of the D° candidate; 

• the cosine of the helicity angle 9* measured between the direction of the charged pion in 
the D° rest frame and the D° direction in the laboratory frame; 



the signed probability Wx as defined in [17| that a given track is compatible with the 



particle hypothesis X, based on its specific energy loss and its measured momentum; 
the signed decay length d xy in the xy plane defined as the distance between the primary 
vertex and the secondary vertex formed by the D° decay products; and 
the largest longitudinal momentum p\ rag relative to the D° flight direction of any track 
which is not used in the D° reconstruction and whose charge is inconsistent with that 
of an accompanying hadron of the D° candidate. This variable helps to separate signal 
candidates in Z° — > cc events from background, since in the formation of D° mesons 
in Z° — > cc events, there is a correlation between the charge of the primary quark and 
that of the fragmentation particle with the highest longitudinal momentum |Tj|. The 



track with the largest longitudinal momentum whose charge is inconsistent with being 
the accompanying hadron is expected to be softer for signal than for background. 

Contributions from D* + decays are suppressed by searching for a track that could have been 
the pion in a D* + — s>D°7r + decay. The corresponding D° candidate is rejected if a track is found 
with the correct charge and a mass difference within 141 MeV < 7711)%+— Wd < 152 MeV. 

A loose preselection of D° candidates is done using the invariant mass m D o, the scaled energy 
xdo, the helicity angle cos 9*, and the particle identification probabilities. The exact cuts are 
listed in table [T]. 

Photons are reconstructed either as showers in the electromagnetic calorimeter, or through 
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conversions 7 



e 1 e 



conversions with the algorithm described in [19 



For the latter, pairs of oppositely charged tracks are identified as 

The z components of the track momenta are 



determined in a fit to the conversion point to improve the momentum resolution. The total 
photon momentum is taken as the sum of both track momenta. 



In the barrel region of the electromagnetic calorimeter, photons are identified as described 
in |Tj| by fitting electromagnetic showers to energy deposits which are not associated with 
any charged particle track. The lateral shower profile is fixed to the Monte Carlo expectation, 
and the normalisation gives the shower energy. The fitted shower energy is then corrected for 
losses in the material in front of the calorimeter using information from the presampler and 
time-of-flight system. Each of the fitted showers is treated as a photon candidate. The photon 
momenta are calculated assuming that the photons originated at the primary vertex of the 
event. 



Any pair of photons (showers or conversions) with an invariant mass between 60 MeV and 
280 MeV is considered as a 7r° candidate. The resolution of the ir° energy is improved by a 
constraint to the nominal ir° mass [[?IJ using a kinematic fit [W \. 



A number of variables similar to those in reference [IS] are calculated for each photon and 
7T° candidate: 



• the energy of the photon or 7r° candidate, denoted E 1 or E n o, respectively; 

• for photon conversions, the signed electron identification probabilities W e calculated from 
the energy loss of the two conversion tracks; 

• for showers in the calorimeter, a variable which parametrises how well the fitted shower 
describes the measured energy distribution in the 9 calorimeter blocks around the fitted 
shower maximum; and three variables that describe their effective separation from the 
closest neighbouring shower and the closest charged track entering the calorimeter. 

The following variables are only used in the selection of 7r° mesons: 

• the invariant mass of the 77 system and the opening angle between the two photons; as 
well as 

• the number of additional photon candidates in a cone around each of the two photon 
candidates under consideration. In each case, the opening angle of the cone is twice the 
angle between the two photons. 

A loose preselection of photons and tt° candidates is done using the photon or tt° energy 
and the invariant 77 mass before the constraint to the n° mass, as listed in table [E]. 

To form D*° candidates, each of the preselected D° candidates is in turn combined with 
every preselected shower photon or 7r° candidate in the event. Figure [I] shows the distribution 
of these D*° candidates in four of the most powerful selection variables. Tighter requirements 
are then applied to further reduce the background level in the D*° sample. They are based on 
the variables listed above and on the total number of D*° candidates in the hemisphere. For 
each variable Xj, a purity function Aj(xj) is calculated based on the simulation as the ratio of 
c^D*° signal candidates to all candidates: 

\ 1 \ si S nal / \ m 

H x i) = t- iii i -jt(zi) , (1) 

(signal + background) 
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D reconstruction 

1.81 GeV < m D o < 1.93 GeV = 
0.3 < x D o 
-0.85 < cos6* < 0.85 
W K (K) < -1% or + 3% < W K (K) 

-35% < W n (K) < 0% 
Wjjr) < -1% or + 1% < W n (ir) 
photon and 7T° reconstruction 

800 MeV < £ 7 or 800 MeV < 
60 MeV < m 77 < 280 MeV 

Table 1: The preselection cuts for the reconstruction of D*° candidates. The selection variables 
are defined in the text. 



and a likelihood is constructed as 

(n \ 1 /n 
Y[Xi{Xi)J . (2) 

Finally, D*° candidates are selected using a cut on the likelihood value. This cut has been 
chosen such that the statistical error on the result for the hadronisation fraction f(c— >D*°) is 
minimised according to an independent sample of simulated events. The distributions of the 
likelihood functions £ are shown in figure ^| together with the cut values. 

In figure ^ the Am=mD.o- m^o mass difference distributions of the selected candidates 
are shown. The peaks at low mass differences are due to the signal from D*° decays. The 
background shape is determined from the data as described below in section [3.2| . Using this 
shape, the D*° yield is extracted from a fit to the mass difference distributions, which is discussed 
in section |4.2| . 



3.2 Background Determination 

Even after applying the likelihood selection, the background level in the samples of D*° candi- 
dates is still high. To minimise the dependence on Monte Carlo modelling, a method has been 
developed to determine the shape of the background in the Am distributions from the data. 

A sample of candidates enriched in background rather than signal is prepared with a method 
similar to the one for the signal selection (see section |3.1|) . The likelihood is modified to select 



background candidates from the m^o sideband regions 1.70 GeV < m^o < 1.84 GeV and 
1.91 GeV < m^o < 2.00 GeV. The sidebands are chosen to lie close to the signal region so that 
the kinematical properties of candidates in the sidebands do not differ significantly from those 
of candidates in the signal region. The modified likelihood function is of the form 

C\x 1 ,...,x n ) = \X' mDO (m D o) X' cosd ,(cos9*) JJ A»(a?i) I , (3) 

\ Xi^m D o, cos 6* ) 

with the A functions left unchanged for all selection variables except m^o and cos 6*. Because 
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the original likelihood function C depends on m D o and cos6**, new functions A^ DO (m D o) and 
^cos6»*( cos ^*) have to be introduced such that the shapes of the distributions of the likelihoods 
C and CJ agree for background. In practise, this is done by constructing the A' functions for 
m D o and cos#* separately, with the same set of values as for the corresponding A function. 
For the selection of background candidates, the same cut is placed on £ as is on C for signal 
candidates. 

One of the main requirements for the background selection is that it correctly reproduces 
the Am shape of true background. If the selection variables m^o and cos 9* are assumed 
to be uncorrelated with the other selection variables and the mass difference Am, the Am 
shapes of background candidates from the same source in the signal and background samples 
agree by construction. For true background, the correlations between m D o (cos^*) and any 
other selection variable have been found in the simulation to be < 1% (< 10%), where the 
largest correlations are between cos#* and p\ rag (—9%) or W n (n) (—7%). It has been tested 
in the simulation that even with this level of correlations, the background shape is correctly 
reproduced for background candidates from the same source. 

Background can be classified into candidates where a correctly reconstructed D° is combined 
with a photon or tt° that does not come from a D*° decay and other background, in which the 
charged pion and kaon candidates do not come from the same D° decay. For these two contri- 
butions, the shapes of the Am distributions are found to be slightly different in the simulation. 
As outlined above, it is expected that the background determination procedure yields the cor- 
rect background shapes individually for both contributions. However, the background sample 
contains fewer candidates with a correctly reconstructed D° than the signal sample. This leads 
to a small bias in the overall background shape derived from the background selection, which is 
taken into account by reweighting this shape according to the sample composition in the signal 
and background samples as found in the simulation. 

Candidates that pass the cuts on C and also on £ are rejected from both the signal and the 
background samples. In the simulation, it has been found that this requirement rejects 8% of 
the true signal from the signal sample, and 27% of the true signal in the background sample. It 
has also been verified that the fraction of candidates removed by this requirement is the same 
within errors for the data and the simulation. 

Although the background sample is depleted in signal candidates relative to the signal 
sample, there is a remaining D*° signal contamination in the background sample. Typically, the 
signal fraction is a factor of 5 smaller (cf. table ^) in the background samples. These candidates 
have the effect of reducing the number of signal candidates after background subtraction and 
are taken into account by calculating an effective efficiency, which is described in the next 
section. 



3.3 Reconstruction Efficiency 

The efficiencies of the D*° reconstruction are determined from a Monte Carlo sample which is 
statistically independent from the one that has been used for the determination of the likelihood. 
In this sample of 6.5 million simulated hadronic events, the numbers of true D*° mesons that 
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pass the cuts on C are computed. The resulting efficiencies in the different decay modes are 
listed in table separately for Z° — > cc and Z° — > bb events. 



As discussed in section |3.2| , some D*° mesons are also contained in the background-enriched 
samples. This number is determined from the same set of simulated events, and an effective 
efficiency is computed from the values given in table |^ that takes into account the relative signal 
to background ratio of the background sample that is given in the last column of the table. In 
the fit to the mass difference distributions described in section |4.2| , these effective efficiencies 
are used. 



decay chain 


distribution 


efficiency, c— >D* 


efficiency, b^D* 


(f) bks /(f) slg 


U K-7T + 


Am(D°,7) 


(2.88 ±0.16)% 


(2.15 ±0.19)% 


0.20 ±0.03 


U K-7T + 


Am(D°,7) 
Am(D°,7r°) 


(1.80 ±0.10)% 
(1.18 ±0.08)% 


(1.01 ±0.10)% 
(1.05 ±0.10)% 


0.20 ±0.03 
0.16 ±0.01 



Table 2: The reconstruction efficiencies for a given decay chain in a given mass difference 
distribution as determined from the simulation. Note that D*° — ► D°7r° decays are measured 
in both the Am(D°,7) and Am(D°,7r°) distributions. The efficiencies are based on candidates 
with scaled energies £d° > 0.3. Only statistical errors are listed. The signal to background 
ratios (f) bkg in the background samples relative to the ratio (f) slg in the corresponding signal 
sample are given in the last column. 



4 Measurement of /(c— ►D* ) 

The dominant sources of D*° mesons in Z° decays are the production in Z° — > cc events and in 
bottom hadron decays. For scaled energies x^a > 0.3, the production of D*° mesons in gluon 
splitting processes g— >cc is highly suppressed and will be neglected in the following. 

The hadronisation fraction /(c— >D*°) is determined in a fit to the Am=mD*o-mDO distri- 
butions of the selected signal and background samples. A simultaneous fit is performed to the 
Am distributions of both D°7 and D°7r° candidates. In the fit, the contributions from b hadron 
decays are subtracted, and the fitted signals are corrected for the efficiency in Z° — > cc events. In 
this section, the determination of the signal contributions from b hadron decays, the fit to the 
mass difference distributions, and systematic uncertainties of the measurement are discussed. 



4.1 Subtraction of the Component from b Hadron Decays 

The determination of the contribution from Z°— >bb events is based on a previous measurement 
of D* + production M. Assuming isospin invariance (i. e., equal production of D*° and D* + 
mesons in b hadron decays, as well as equal hadronisation fractions /(b— >B~) = /(b— >-B ) in 
Z° — >bb decays), the measured D* + production rate in Z° — >bb events is used as an estimate 
of the corresponding D*° production rate. This is done for scaled energies x^o > 0.3, and the 
remaining c— >D*° signal is then extrapolated to the full xd° range. 
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Previously, OPAL has measured the production of D* + mesons in Z° — > bb events to be 
R h f (b^D*+) Br (D* + ^D°tt + ) Br (D°^K-n + ) = (1.334 ±0.049 ±0.078) x 10" 3 §, where R h 
is the partial hadronic decay width of the Z° boson into bb pairs. Using this number and the 
branching fraction Br(D* + ^D°7r+) = 0.683 ± 0.014 [|U], the D*+ production rate in Z°^bb 
events for x D o>0.3 is calculated to be 

/3 = R h f (h^D* + ) Br (D°^K-7i + )\ , = (0.980 ± 0.068) x 10~ 3 . (4) 



The error includes both the statistical and systematic errors. The restriction to £d°>0.3 is based 
on Monte Carlo simulation and the measured parameters of the b fragmentation function. 

This information is used to obtain the number of reconstructed D*° signal candidates from 
b hadron decays, which is then subtracted from the fitted signal. Thus, the hadronisation 
fraction f(c— >D*°)\ „ „ is determined in the fit from 

J v ' lx D o>0.3 



x D o>0.3 



iy sig (D^D°X) - (3Br(D*°^D°X) £ ^° x 

2iV had /_x 

D*°->D°X ' \°' 

X=7,7T° C 

where R c is defined in analogy to R^, N sig denotes the signal as obtained from the mass difference 
distribution, iVhad is the number of hadronic Z° decays analysed, e c and are the effective 
efficiencies for reconstructing D*° decays with x D o>0.3 from Z° — > cc and Z° — > bb decays, 
respectively, and the sum runs over the two decay modes of the D*°. 

The assumption of equal D*° and D* + production in b hadron decays is justified from a 
CLEO measurement ^TJ of 

) J - = 1.03 ±0.14. 6 

5r(B^D*+X) V ; 

The error on this measurement is taken into account as described in section |4.4.4| . Decays 



of B and B~ mesons account for most D* mesons from b hadron decays at LEP, since B s 



mesons preferentially decay to final states with a [ 20j , and b baryons are expected to decay 
dominantly to final states with a c baryon. 



4.2 The Fit to the Mass Difference Distributions 

The number of reconstructed D*° decays is determined in a simultaneous fit to the mass dif- 
ference distributions of the selected signal and background samples in the two decay modes 
D*°^D°7 and D*°^D°7r . 

Two signal peaks are visible in the mass difference distribution of D°7 candidates, which is 
shown in figure ^a. The peak at the nominal D*°— D° mass difference of 142 MeV p(J is due 
to D*° — > D°7 decays, whereas the one at lower Am values is from D*° — > D°7r° decays, where 
only one photon from the ir° decay is combined with the D° candidate. 

In the mass difference distribution of D°7r° candidates, the signal from D*° — > D°7r° decays 
is expected around the nominal mass difference. Additional contributions are expected from 
D*° — > D°7 and D*° — ► D°7r° decays, where only one of the two photons from the 7r° decay is 
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used. In both cases, an unrelated photon candidate is added to form a n° candidate. In the 
simulation, both contributions have been found to be broad and very similar in shape to the 
overall background. Therefore, they are not explicitly accounted for. 

In both mass distributions, the signal is parametrised with a functional form, while the 



background shape is determined from the background sample as decribed in section |3.2| , with 
its normalisation determined in the fit. A sum of two Gaussians is used as signal parametrisation 
for the Am distribution of D°7 candidates (see figure |3|c and d) in order to accomodate the 
contributions from both D*° — > D°7 and D*° — > D°7r° decays. The mean values of the Gaussians 
are fixed to 145 MeV (D*° — > D°7) and 80 MeV (D*° — > D°7r°) as determined from the simulation, 
and the widths are allowed to vary. For D°7r° candidates, a modified Gaussian of the form 

dn f \ ( 1 /Am — m 7T o\ 2 \ 

Am — m n o J exp — - I I ; Am > m^o (7) 



d(Am) ^ * J 1 \ 2 V o 

is taken as the signal parametrisation (see figure ||g and h) to account for threshold effects. 
Here, m n o denotes the nominal tt° mass f2"0fl , and the parameter a and the normalisation are 
determined in the fit. In the simulation, such a function has been found to parametrise the 
distribution of D°7r° signal candidates well. 

The hadronisation fraction f(c — > D*°) is determined from a simultaneous fit to both Am 
distributions, with the constraint that the efficiency corrected number of D*° — > D°7r° de- 
cays be the same whether determined from the peak at lower Am values in the Am(D°,7) 
distribution or from the Am(D°, tt°) distribution. Also, the efficiency corrected numbers 
for the two D*° decay channels are fixed to the world average value of the branching ratio 
Sr(D* Q ^D 7)/Sr(D*°^D Q 7r ) = 0.616 ± 0.076 Q, and it is assumed that these two decay 
modes saturate the D*° width. The contribution from bottom hadron decays, which is derived 



from an earlier OPAL measurement as described in section 4.1, is subtracted in the fit. 



The fit results are illustrated in figure |3|, where the obtained signal parametrisations are 
shown together with the background subtracted mass difference distributions. The combined 
fit for the D*° measurement has a \ 2 °f 61.6 for 69 degrees of freedom. The hadronisation 
fraction /(c — > D* )| x . o>03 is extrapolated to the full range of scaled energies x^o using the 
simulation; it is found to be 

R c f(c^D*°)Br(D°^K-iT + ) = (1.44 ± 0.36) x 10" 3 , (8) 
where the error is statistical. 



4.3 Consistency Checks 

A number of consistency checks has been performed, in particular to test the sensitivity of the 
result to the background subtraction procedure. Monte Carlo simulation has been used to check 
that the procedure as outlined in section [T2] accurately describes the shape of the background 
in the signal sample. Also, it has been verified that the shape is not sensitive to the specific 
choice of A' functions. The background determination is sensitive to correlations between moo 
or cos 6* and any other selection variable. Any correlation seen in the data is well reproduced 
in the simulation. Finally, the analysis has been repeated on the simulation, and the generated 
rates are reproduced within the statistical errors. 
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For all selections, the likelihood cut has been varied and the analysis repeated. The ranges 
of cut values correspond to relative changes in efficiency of about ±50%. The variations have 
been found to be consistent with statistical fluctuations. 

In the fit, the ratio Br(D*° — > D°7) / Br(D*° — ► D°7r°) is fixed to its world average value 
of 0.616 ± 0.076 |2g. If instead, it is allowed to vary, a value of Br(D*° -> D°j)/Br(D*° -»■ 
D 7T°) = 0.39 ± 0.17 (stat.) is found, which is consistent with the world average. The resulting 
hadronisation fraction /(c^D*°) changes by —7.5%. 



4.4 Systematic Uncertainties 

Systematic uncertainties are related to the modelling of the selection variables, the detector 
resolution, the procedure used to subtract the background, the determination of the b — > D*° 
contribution to the measured signal, and to the extrapolation to the full range of scaled energies 
xdo. In the following sections, each of these categories is discussed in turn. The relative values 
of all errors are listed in table 



4.4.1 Uncertainties from the Modelling of Selection Variables 

Possible differences between distributions in data and simulation could influence the efficiencies 
and the background determination. These effects are studied separately for each selection 
variable. The resulting systematic errors are assumed to be uncorrelated and are therefore 
added in quadrature. 

Two principal methods are used to determine systematic uncertainties in the modelling of 
a selection variable: 

(A) Variables for which the distributions of signal candidates are measured are treated by 
reweighting the events in the Monte Carlo simulation such that for the weighted events, the 
simulated signal distribution agrees with the measured one. Using the likelihood given in 
equation (H), the background determination and the fit are then repeated and the original result 
is corrected according to the observed difference; the uncertainty in the difference is interpreted 
as a systematic error. This method takes correlations between the selection variables into 
account to the extent that in the case of non-zero correlations, the reweighting of events alters 
also the distribution of any variable that is correlated with the one variable in question. This 
procedure is applied 

• to the invariant mass distribution of D° candidates, 

• to the Peterson et al. fragmentation parameters which have been varied in the ranges 
corresponding to mean scaled energies of D*° mesons in Z° — > cc decays of 0.506 < 

zo-^ C c^D* < 0.531 and of weakly decaying b hadrons in Z° — ► bb decays of 
0.702 < (x(X b )) z o^ b b^x b < 0.730 |22|, respectively, 



to the lifetime of D° mesons, which has been varied within rpo = (0.415 ± 0.004) ps f20|j . 



to the distributions of energies E 1 and E n o, where the simulated events are reweighted 
such that the helicity angle distributions in D*° — > D°7 and D*° — > D°7r° decays agree with 
the measured D* + spin alignment in Z° — > cc events (the spin density matrix element is 
poo = 0.40 ±0.02 pi), 
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relative statistical error 


24.9% 


relative systematic errors: 




modelling of selection variables: 




D° mass resolution 


1.3% 


(x(D*)) z o_ +cg ^ D »o 


3.9% 


(x(X b )) z o^ bS ^ Xb 


3.0% 


D° lifetime 


0.8% 


D*° spin alignment 


0.3% 


effective shower isolation 


9.7% 


photon pair opening angle 


2.8% 


number of n° candidates 


0.2% 


number of D*° candidates 


1.3% 


dE/dx 


3.6% 


dE/dx preselection cuts 


3.1% 


fragmentation tracks 


2.7% 


shower fit 


5.9% 


detector resolution: 




tracking resolution 


6.4% 


calorimeter energy scale and resolution 


1.7% 


fit procedure: 




background normalisation 


3.4% 


Br (D*° -> D° 7 ) /Br (D*° -> D tt ) 


2.0% 


background shape 


5.5% 


contributions from other D° decays 


0.4% 


contributions from D ^ and Dg decays 


1.4% 


b subtraction and extrapolation: 
subtraction of the b — > D*° contribution 
extrapolation to x^o = 


12.2% 
1.2% 


total relative systematic error 


20.9% 



Table 3: A breakdown of the relative statistical and systematic errors on the hadronisation 
fraction /(c^D*°). 

• to the effective separation of a shower from the closest other shower in the electromagnetic 
calorimeter, 

• to the opening angle of a pair of photons, and 

• to the numbers of reconstructed n° and D*° candidates. 

The last three of the above variables depend on the event topology. For these variables, possible 
deviations between data and simulation affect the signal and the background in the same way. 
The simulated events are therefore reweighted such that the distributions of all candidates in 
data and simulation agree. The analysis is then repeated, and the observed difference in the fit 
result is treated as a systematic error. 

(B) The second method is used for the remaining selection variables, which are 

• the track probabilities calculated from dE/dx information, 

• the largest longitudinal momentum of fragmentation tracks inconsistent with being an 
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accompanying hadron, and 
• the goodness of the shower fit. 

The corresponding function Aj(xj) is set to its maximum value for one variable Xj at a time. The 
resulting modified likelihood functions Li and L\ are independent of that particular variable 
Xi, so using these functions instead of the original likelihoods removes any possible bias due 
to the simulation of x%. The selection of the signal and background D*° candidate samples is 
repeated with these modified likelihoods. By construction, these samples also contain those 
candidates selected with the original likelihoods £ and £'. The ratio of the fit results obtained 
with the original (C,C) and the modified likelihoods is calculated for the data and 

the simulation, where the statistical correlation between the samples is taken into account. 
The relative difference of these ratios is interpreted as the systematic error associated with 
the variable In the case of the signed particle identification probabilities, the cuts for the 
preselection of charged kaons and pions (see section |3.1| ) are retained in order to obtain a clear 
signal. In reference [|[, an error of 3.1% has been determined for these cuts in the D° selection. 
This error is included as an additional systematic uncertainty. 



4.4.2 Uncertainties in the Detector Resolution 



Tracking resolution: 

Uncertainties in the modelling of the central detector are assessed by repeating the analysis 
with the tracking resolutions varied by ±10% around the values that describe the data best. 
The redetermined efficiencies are compared with the original ones, and the relative difference 
is interpreted SIS db systematic error. 

Calorimeter energy scale and resolution: 

The energy scale and resolution of the electromagnetic calorimeter and the multiplicity and 
resolution in the time-of-flight scintillators are treated analogously. The reconstructed 7r° mass 
distribution has been used to determine the corresponding resolution and scaling parameters. 
Thus, width and position of the n° mass peak are well reproduced in the simulation |19| , and 
any possible bias is accomodated by the variation of the detector resolution. 



4.4.3 Uncertainties in the Fit Procedure 



Background normalisation: 

To check the determination of the background normalisation, the fit to the mass difference 
distributions has been repeated with the range restricted to values of Am(D°,7)<0.4 GeV 
and Am(D°, 7r°)<0.2 GeV. Deviations from the previous results were interpreted as systematic 
errors. 

D*° branching ratio: 

The fit to the Am(D°,7) and Am(D°,7r°) distributions is constrained to the world average of 
the branching ratio Br(D*° — > D°7)/_Br(D*° — > D°7r°) = 0.616 ± 0.076 @. The fit is repeated 



with this ratio varied within its errors, and the observed difference is taken as a systematic 
error. 
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Background shape: 

Background candidates with a correctly reconstructed D° lead to a bias in the Am distribution, 



which is taken into account in the background subtraction as outlined in section |3.2| . The 
analysis has been repeated without correcting for this bias, and half the difference from the 
original result is assigned systematic error. 

The origin of background candidates has been studied in the simulation, and no other source 
for a bias of the background shape has been identified. In particular, it has been checked that 
the contribution from photons or tt° mesons from D**, B*, and B** decays is small and does 
not exhibit a pronounced structure. 

Contributions from other D° decays: 

In the simulation, the contribution to the signal sample of D*° candidates from D° mesons not 
decaying to K~n + was found to be 5.9%. This contribution is accounted for in the efficiency 
determination. However, mismodelling of the branching fractions for these other decay modes 
can introduce a potential bias in the D*° measurement. The corresponding systematic error is 
evaluated assuming relative contributions of other decay modes as in |§ and assigning errors 



according to the errors on the branching fractions as given in [ ZD] . 



Contributions from D*+ and D*+ decays: 

In the reconstruction of D*° decays, there are contributions from D* + and D* + mesons, which 
can also decay via 7 or tt° emission. In principle, these decays can lead to a signal in the Am 
distributions similar to that to be measured. However, it has been found in the simulation 
that the relative contribution from D* + and D* + decays to the selected signal and background 
samples is approximately equal. A small correction is applied to account for residual effects, 
and half this correction is assigned as a systematic error. 

It has been checked that the statistical correlation between the D*° — > D°7r° signals in the 
Am(D°, 7T°) and Am(D°, 7) distributions is negligible. Also, there is only a negligible fraction 
of cases where both photons from a D*° — > D°7r° decay lead to an entry in the Am(D°,7) 
distribution. 



4.4.4 Uncertainties from the b Subtraction and the Extrapolation 



Subtraction of the b— >D*° contribution: 

As mentioned in section 4~T, it is assumed that the production of D*° and D* + mesons in 
b hadron decays is equal for scaled energies xd°>0.3. The error on the production of D* + 
mesons in Z° — > bb decays as given in equation (|J) is taken into account. In addition, the 
production of D mesons in B and B meson decays is varied within the range of the CLEO 
measurement given in equation (|6]) to assess the systematic error due to this assumption. At 
LEP, W and B" mesons account for (75.6 ± 4.4)% of all weakly decaying b hadrons [pOfl . To 
evaluate the production of D*° mesons in decays of other bottom hadrons, the world average of 
Br(B° or B ^D*+X) = (23.1 ±3.3)% @ is compared with the recent OPAL measurement of 
J R b /(b^D*+) J Br(D* + ^D 7r + ) J Br(D°^K-7r + ) = (1.334 ± 0.049 ± 0.078) x 10~ 3 §. From the 
branching fractions as given in |2(| , the production of D* + mesons in B^ and A b decays is then 
found to be Br(B>l or A b ^D* + X) = 0.25 ± 0.13. Assuming equal D*° and D* + production in 
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B" and A b decays, this leads to R h f(b^>W B or A b -^D*°X) = 0.013. A 100% error is assigned 
to this quantity. 

Extrapolation to x D o = 0: 

The extrapolation of the measured quantity R c /(c— >D*°) Br(D° — >K~7r + ) \ x >Q 3 from x D o>0.3 

to the full range of scaled energies is based on the fragmentation function of Peterson et al. [15[ . 
To assess the uncertainty associated with this extrapolation, it is repeated with mean scaled 
energies of D*° mesons from Z° — > cc decays varied in the range given in section 4.4.1 . The 
difference from the previous result is interpreted as a systematic error. 



4.5 Results of the D*° Measurement 

The production of D*° mesons in Z° — > cc events is measured to be 

J R c /(c^D*°) J Br(D ^K-7r+) = (1.44 ± 0.36 (stat.) ± 0.30 (syst.)) x 10~ 3 . (9) 

From this value, from the branching fraction Br(D° — >K~n + ) = (3.83 ± 0.12)% ||20|| , and from 
the standard model prediction of R c = 0.172 p^| , the hadronisation fraction f(c — > D*°) is 
computed to be 

/(c^D*°) = 0.218 ± 0.054 (stat.) ± 0.045 (syst.) ± 0.007 (ext.). (10) 
Here, the last error corresponds to the error on the D° branching fraction. 



5 Measurement of D* + Production 



The decays of D* + mesons are very simliar to those of D*° mesons. The techniques developed 
for the D*° analysis can therefore be used for the reconstruction of D* + mesons. 

The D* + meson decays dominantly to the D+7 final state, since the D+7r° channel is sup- 
pressed by isospin invariance: Br(D*+ -> D+n )/ Br(D* s + -> D+7) = 0.062±g;g?g ± 0.022 E5 



Thus, for the measurement of the production rate, only the D* + — > D+7 decay is used. The 
D + mesons are reconstructed in their decay chain — > 07r + , — > K + K _ . Candidates with 
invariant masses m^,<1.05 GeV and 1.90 GeV<m D +<2.04 GeV, a scaled energy of x D +>0.35, 
and a helicity angle of the pion in the rest frame satisfying — 0.90<cos#*<0.95 are retained 



if the dE/dx requirements stated in section |3.1| are fulfilled for the kaon and pion candidate 



tracks. The preselection of photons is identical to that described in section 3.1 



The D* + candidates are selected with a cut on a likelihood using the variables described in 



section [3.1| accordingly, and taking the reconstructed mass of the <p and the cosine of the angle 
between the and one of the kaons in the 4> rest frame as additional inputs. The background 
shape is determined from candidates in sidebands of and m D + using the technique described 



in section |3.2|. The resulting mass difference distribution is shown in figure 



No attempt is made to separate the contributions from the processes c— >D* + and b— »D* + . 
However, the efficiencies for D* + reconstruction in Z°— s>cc and Z°^bb events have been found 
to be equal within errors. Thus, it is still reasonable to extract the overall D* + production in 
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Z° decays from a fit to the Am(D+,7) distribution. No isospin violating decays D* + ^D+7r° 
have been generated in the simulation; in the data, a small contribution from such decays is 
expected at Am values below the nominal D* + -D+ mass difference as shown in figure |]c. A 
double Gaussian is used to parametrise the signal, where the Gaussian at lower Am is fixed 
to the expectation and varied by ±100% in case of the data. The number of D* + mesons per 
hadronic Z° decay with x D +>0.35 is found to be 

n(Z°^D:+)Sr(D*+^D+ 7 )Sr(D+^07r+)Sr(0^K+K-)| a;D+>Q35 

= (7.1 ±L9 ±1.7) x 10" 4 , 

where the first error is statistical and the second systematic. The fit to the data has a x 2 °f 
37.1 for 43 degrees of freedom. 

The contributions to the systematic error were evaluated in a similar way as described 
earlier. Table ^ lists the most important errors. The largest systematic error on the value in 
equation (|i~I[) is introduced from the uncertainty in the background shape. The background 
sample is taken from a sideband in m^, and a correction has to be applied to account for the 
fact that fewer correctly reconstructed <f)—>K + K~ decays enter the background than the signal 
sample. When the fit is repeated with the same signal parametrisation plus an exponential to 
account for potential problems in the modelling of this bias, a consistent result is found. The 
deviation from the previous result is interpreted as systematic error. No other sources have 
been identified that could lead to a significant bias in the background shape. 



relative statistical error 


27.4% 


relative systematic errors: 




extrapolation to x D +=0 


(19.6%) 


background shape 


11.8% 


tracking resolution 


8.5% 


D+ lifetime 


8.5% 


effective shower isolation 


8.1% 


shower fit 


6.9% 


dE/dx preselection cuts 


6.4% 


(z(Xb)) Z o_ > bb-»X b 

others (evaluation similar to the D*°) 


6.2% 


9.5% 


total relative systematic error 


23.8% 
(30.8%) 



Table 4: A breakdown of the relative statistical and dominant systematic errors on the D* + 
measurement. The errors in brackets apply only when considering the value which has been 
extrapolated to all scaled energies x B +. 



The above D* + measurement is extrapolated to the full range of scaled energies x D + . Here, a 
large systematic uncertainty is introduced from the unknown relative contribution from c— >D* + 
and b — > D* + components. When varying the unknown ratio f(b — > D* + )//(b — > D+) within 
0.6 ± 0.2, one finds the total rate 

n(Z°^D*+) 5r(D:+^D s +7) Br(D+ - 0tt+) flr(0-K+K-) (m 

= (1.69 ± 0.46 ± 0.40 ± 0.33) x 10~ 3 , 1 ' 
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where the third error is introduced from the extrapolation to the full range of scaled energies. 
This rate is consistent with expectations and can be regarded as a further cross-check of the 
D*° analysis, which uses very similar techniques. 



6 The Relative Production Rate of Charmed Vector 
Mesons in Z° — ► cc Decays 



The relative production rate of vector mesons containing the primary quark, Py, is an important 
parameter in fragmentation and hadronisation models. Hadrons that contain the primary quark 
can in principle be produced either directly or in decays of higher resonances. Therefore, there 
are two possible definitions of Py: the quantity P v s , where inclusive production (including 
decays of higher resonances) is considered, and Py im , which is defined when considering direct 
production only. These two values may differ due to the effects of decays of higher resonances, 
where some decays are forbidden by spin and parity conservation. 

In previous investigations of Py in the charm system 0, f|, ||, isospin invariance between 
the neutral and charged non-strange vector mesons was assumed. Given the measurement of 
the hadronisation fraction f(c— >D*°) presented in the first part of this paper, an explicit check 
of this assumption is now possible. 

Recently, a measurement of the spin alignment of D* + mesons in Z° — > cc decays has been 
presented [p3[ . While Py denotes the relative production of vector and pseudoscalar mesons, 



spin alignment measurements provide information on the relative production of different vector 
meson spin states. Thus, the combination of these two measurements provides further insight 
into the inclusive production of charmed mesons in the fragmentation process j3|. 

In the first part of this section, the existing OPAL measurements of the production of 
charmed mesons with no orbital angular momentum are used to derive a value of P{f, which 
is interpreted in conjunction with the D* + spin alignment. 

In the second part, additional input from excited D states is used to derive Py im . This 
quantity can be compared more directly to model calculations, since such models generally make 
predictions for primary hadron production. However, a determination of Py im is experimentally 
challenging because of the difficulties in assessing the fraction of hadrons produced in decays 
of higher resonances. In this paper, the production of L=l mesons is taken into account, while 
production and decay of higher states is considered a part of the fragmentation process. For 
each of the light flavours u, d, and s, four charmed mesons with orbital angular momentum L=l 
are predicted. The naming convention is Dj , where J denotes the total spin of the meson, and 
an asterisk indicates that the meson has parity (— 1) J . In the following, these are collectively 
referred to as D** mesons. Using the recent OPAL measurements of the hadronisation fractions 
/(c— >D5), /(c — ^Dg ), and /(c— ►D^) fllPl , decays of D** mesons can be taken into account 



to study the pP nm ratio. First, the dependence of Py im on the production of the unmeasured 
L=l resonances is discussed, and a model independent formula for the calculation of Py im is 
derived. Second, the production of the unmeasured L=l resonances is assessed in a simple spin 
counting picture. 
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The following discussions are based on the /(c— >• D*°) measurement described in the first 
part of this paper and on previous OPAL measurements of charmed meson hadronisation frac- 
tions. Table |5| contains an overview of the values entering the following computations. 



charmed meson production, OPAL measurements 



value and reference 



Rc /(C- 

Rc /(c- 
Rc /(c- 
Rc /(c- 



•D*°) Br(D°- 
•D*+) Sr(D H 
■D°) 5r(D°- 
■D+) fir(D+- 



K-7T + ) 



5r(D° — > K~7r" 



7T 



or D*°) 5r(D° or D 

/(c— D°— D*+ti-) 



*0 
2 



D*~ 



7T 



/(c 
i?c/(c 



>D* + 7T-) + /(c- 



+u 2 



>D* + 7T-) 



(1.44 ± 0.36 ± 0.30) xl0~ 3 (section^ 
(1.041 ± 0.020 ± 0.040) x 10- 3 @ 
(0.389 ± 0.027j#{g!) x 10~ 2 § 
(0.358 ± 0.046lHi) x 10~ 2 

0.5+0.2^ 
0.7-0.3, 



(2.8±g;f ± 0.3 ± 0.4) x 10~ 3 [| 



(4.2 d= l.l±g:f±8;^ x io- 3 
0.56 ± 0.1518$ 



branching fractions and ratios 



value and reference 



Br(D + -> K-TT+7T+) /Sr(D° -> K-tt+) 
J Br(D* + ^D°7r+) 



Br(DS-+D*?r) 



Br(D*^D*7r)+Br(D*^D7r) 

_Br(D*+^D*K) 
Br (D*^— >D*K)+Br (D* 9 + ->DK) 



Br(D*^D* 



it 



flr(D£-D*K) 



0.0383 ±0.0012 [gO 
2.35 ±0.16 ±0.16 
0.683 ±0.014 [|0 

0.311 ±0.051 [|0| 
0.107 ±0.016 
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assumption on relative D** meson production in a spin counting model 



Mil 



/(c^D*°) : /(c->D?) : /(c-D?) : /(c-D 
/(c^D: + ):/(c-,D^):/(c-,D s + 1 ):/(c^D: 2 + 



1:3:3:5 



Table 5: The measurements and assumptions entering the computation of values. 



6.1 Tests of Isospin Invariance 

Isospin invariance suggests equal primary production rates for corresponding cu and cd mesons. 
Thus, the hadronisation fractions /(c^D*°) and /(c— >D* + ) are expected to be the same, as 
long as decays of higher resonances contribute equally to both D*° and D* + production. For 
the vector mesons D*° and D* + , isospin invariance can therefore be tested directly from the D*° 
and D* + production rate measurements and the branching fraction Br(D* + — > D°7r + ) as listed 
in table [|, resulting in a ratio of 

/(c^D*°) 
/(c^D*+ 

consistent with 1. 



0.94 ±0.31, (13) 



Since D*° — > D+rr - decays are kinematically forbidden while D* + — > D°7r + decays are not, 
the observed hadronisation fractions of D° and D + mesons are expected to differ even if the 
primary production rates are equal. When taking into account these D* + decays and assuming 
isospin invariance between D° and D + mesons, the ratio 

R = /MH^ (14) 

/(c^D*°) + (2fir(D*+^D°Tr + ) - l) /(c^D* + ) 



is expected to be equal to one. From the values given in table §, one obtains R = 1.19 ± 0.36, 
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consistent with 1. Under the assumption of equal hadronisation fractions f{c— >D*°) = /(c— ► 
D* + ), other experiments Q have obtained values consistent with this result. 

With isospin invariance thus confirmed within experimental errors, the hadronisation frac- 
tions of D*° and D* + mesons are assumed to be equal in the following. The mean D*°/ + 
hadronisation fraction is then defined as the weighted average of the values of f(c — >D*°/ + ) de- 
termined from equation QUI) and from the direct measurements of /(c— >D*°) and /(c— >D* + ). 
The resulting value of 

R c /(c^D* 0/+ ) 5r(D* + ^DV) Pr(D°^Rr7r + ) = (1.045 ± 0.044) x 10~ 3 (15) 

is used in the computation of Py in the following sections. 



6.2 Effective Charmed Meson Production 



The effective value P v s is calculated from the mean D* 0,/+ hadronisation fraction derived in the 
previous section and the measurements listed in table || to be 

pes _ 2/(c^D* 0/+ ) _ n p- 7 + nn p- / lfi N 
F v - /(c ^ D o) + /(c^D+)" U -^ ±U - U5 - {Vb) 
Here, the assumption has been made that isospin invariance is valid, as tested in the previous 
section. If, instead, the measured hadronisation fractions are directly combined without the as- 
sumption of isospin invariance for D*° and D* + production, a value of P y ff = ^fe^o)^jfe§+) — 
0.55 ± 0.10 is obtained. In principle, the measurements are expected to be correlated. The 
largest correlation is expected between f(c— > D°) and f{c— >D* + ) and has been estimated to 
be smaller than 30%. This introduces a systematic error on Pf of less than 0.01. 



The result from equation fllq) can be interpreted in connection with the D* spin alignment. 
In the fragmentation process, four different spin states of charmed mesons without orbital 
excitation can be formed: The vector mesons D*°, D*+, and D*+ with states J = 1 and 
m = —1,0,-1-1, as well as the pseudoscalar mesons D°, D + , and D+ with J = m = 0. The 
relative inclusive production probabilities P for these states are related to Py ff and to the spin 
density matrix element poo of D* mesons from Z°— >cc events via 0, |29 



Pf = V^ + V&- and 

poo = vy^/ivf^ + v?^ 1 ) . 1 } 

Here, P^Jf 1 denotes the sum of the production probabilities for the m = +1 and m = —1 
states that cannot be distinguished experimentally. The spin density matrix element poo gives 
the probability to find a vector meson in the m = state. A simple spin counting model 
suggests values of P^Jf 1 = \ and P^jf = Py= ° = \ (see for instance ||). 



2 ' J=l ' J=0 4 

For D* + mesons from Z°-^cc decays with scaled energies xd*+>0.2, p o has been measured 
at OPAL to be p o — 0.40 ± 0.02 [23]. From this result and the above value of Py ff , the 



production probabilities 

p™ = =° = Poo pf = 0.23 ±0.02, 

p^jf 1 = (1 - Poo )Pf = 0.34 ± 0.03 , and (18) 

Vy= ° = 1-Pf = 0.43 ± 0.05 
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can be derived, where it has been assumed that the measured spin density matrix element poo 
applies to both D*° and D* + mesons in Z° — > cc events. 

These probabilities show a clear deviation from the simple spin counting picture. The 
production of vector mesons is suppressed in favour of pseudoscalar mesons. This is mostly 
due to a suppression of the m=±l vector states, while the production of the J—l, m=0 state 
agrees within errors with a spin counting picture. 

An overall suppression of vector meson production is, for instance, expected in thermody- 
namic models, where the constituents of hadrons are pictured as a gas with a temperature T, 
such that a relative suppression of the heavier states by a factor of exp(— Am/T) is predicted. 
From the above value of P v s , the temperature is calculated to be T = (174 +36) MeV, con- 



sistent with the value determined in [3D|. It should be noted, however, that thermodynamic 



models fail to explain the observed non-zero spin alignment. 



6.3 Primary Charmed Meson Production 

For the determination of Py nm , the effects of the decays of L=l charmed mesons have to be 
taken into account. The hadronisation fractions /(c— ►D"), /(c^D^ ), and /(c— >D^) have 
been measured at LEP, whereas the production of the broad resonances Dq and D[ has not 
yet been measured. Since the charged D** + mesons have not yet been observed in Z° decays, 
isospin invariance is assumed to be valid in charm fragmentation to assess the production of 
these resonances, which implies equal hadronisation fractions for corresponding D** + and D**° 
states. Furthermore, it is assumed in the following that for D** mesons, the relative production 
of the different spin states does not depend on the light quark flavour. 

The ratio P v " m can be expressed as a function of two unmeasured hadronisation fractions, 
/(c— >Dq°) and f(c-^Df), and two known parameters A and B which depend on the measure- 
ments listed in table ||: 

B-R c [/( c -D5°)+/(c-Df 
The complete derivation and the exact formulae for A and B are given in appendix A. The 
parameters A and B are found to be 

A = (2.29 ±0.34) x 1(T 2 and 
B = (3.82 ± 0.89) x 10 

with a positive correlation between A and B of 77%. 

In figure |5], the results for P v Iim are shown as a function of the two unknown hadronisation 
fractions. In general, low values of P v vim are obtained for small Dq and large D[ production and 
vice versa. This analysis shows that the value of P v nm is not very sensitive to the production of 
the unmeasured broad D** resonances. Generally, the range of P v mm values is consistent with 
predictions. However, to test a given model, a clearer statement can be made when following 
the procedure outlined below. 

In contrast to the above discussion, a test of any specific fragmentation model can be 
performed when using its prediction for the relative primary production of the different L=0 



pprim ^ J-^cJ \y *->\ ) /^qx 



(20) 
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and L=l states. While the prediction for the L=0 states is directly equivalent to a prediction 
of Py im , the model prediction in the L=l sector provides information on the production of 
the unmeasured L=l resonances, which can in turn be used in conjunction with the measured 
hadronisation fractions to obtain an experimental value of Py im . Such a test is described in 
the following for the simplest case, a spin counting model. 

A spin counting model predicts the relative primary production of D and D* mesons to be 
1:3. At the same time, for D** mesons, the relative production of the Dq : D[ : Di : D£ 
resonances is predicted to be 1 : 3 : 3 : 5. When using the latter prediction, P v rim is found to 
be 

P£ rim (spin counting for D**) = 0.55 ± 0.08 . (21) 

However, since the validity of spin counting was assumed in the calculation, the value of P v nm 
is fixed and should be 0.75. Thus, from the discrepancy between the above value and 0.75, the 
simultanous description of both L=0 and L=l charmed meson production in a spin counting 
picture is disfavoured by 2.7 standard deviations. Within the framework of this model test, 
the relative production of the four D** spin states need not be taken from experiment, since 
it is predicted by the model. Thus, the experimental error on the model test is reduced as 
compared to the model independent values shown in figure [5]. As mentioned above, the value 
in equation (El]) is part of a consistency check of a specific model and should not be interpreted 
stand-alone measurement. 

In principle, a comparison of the hadronisation fractions /(c^D* + ) and f{c— >D+) yields 
a model independent measurement of P v nm . However, the result from section [5] cannot be 
interpreted in terms of vector and pseudoscalar meson production in charm hadronisation, 
since no information on the separation of the c— >-D* + and b^D* + components exists. 



7 Summary and Conclusion 

A first measurement of the hadronisation fraction /(c^D*°) in Z°-^cc decays is presented: 

/(c^D*°) = 0.218 ± 0.054 (stat.) ± 0.045 (syst.) ± 0.007 (ext.) . 

This result is consistent with the expectation from isospin invariance. The production rate of 
D* + mesons in hadronic Z° decays has been measured for the first time: 

n(Z°->D*+) Br(D* + ^D s + 7 ) 5r(D+->07r+) £r(</>^K + KT) 

= (1.69 ± 0.46 (stat.) ± 0.52 (syst.)) x 10~ 3 . 

The relative production of vector charmed mesons in Z° — > cc events, Py, is evaluated 
both considering inclusive production (P v s ) and taking into account the effects of secondary 
production in D** decays (Py im ). A value of 

Pf = 0.57 ±0.05 

has been derived from OPAL measurements, consistent with previous results. The dependence 
of Py im on the unmeasured D** multiplicities is determined in a model independent calculation, 
where a weak dependence on the production of the unmeasured broad D** resonances has been 
found. From the determination of P v nm , it is found that for the description of the production 
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of charmed mesons in Z° — > cc decays, a simple spin counting picture is disfavoured by 2.7 
standard deviations. 



Acknowledgements: 



We particularly wish to thank the SL Division for the efficient operation of the LEP accelerator 

at all energies and for their continuing close cooperation with our experimental group. We 

thank our colleagues from CEA, DAPNIA/SPP, CE-Saclay for their efforts over the years on 

the time-of-flight and trigger systems which we continue to use. In addition to the support staff 

at our own institutions we are pleased to acknowledge the 

Department of Energy, USA, 

National Science Foundation, USA, 

Particle Physics and Astronomy Research Council, UK, 

Natural Sciences and Engineering Research Council, Canada, 

Israel Science Foundation, administered by the Israel Academy of Science and Humanities, 

Minerva Gesellschaft, 

Benoziyo Center for High Energy Physics, 

Japanese Ministry of Education, Science and Culture (the Monbusho) and a grant under the 
Monbusho International Science Research Program, 
German Israeli Bi-national Science Foundation (GIF), 

Bundesministerium fur Bildung, Wissenschaft, Forschung und Techno logie, Germany, 
National Research Council of Canada, 
Research Corporation, USA, 

Hungarian Foundation for Scientific Research, OTKA T-016660, T023793 and OTKA F-023259. 



24 



References 

[1] K. Cheung, Recent Progress on Perturbative QCD Fragmentation Functions, Proceedings 
of the 1995 PASCOS/Hopkins Workshop, Baltimore, or hep-ph/9505365; 
E. Braaten, K. Cheung, S. Fleming, and T.C. Yuan, Phys. Rev. D51 (1995) 4819; 
Y.-Q. Chen, Phys. Rev. D48 (1993) 5181; 
M. Suzuki, Phys. Rev. D33 (1986) 676; 
R. Suaya, J.S. Townsend, Phys. Rev. D19 (1979) 1414. 

[2] OPAL Collaboration, K. Ackerstaff et al, Z. Phys. C74 (1997) 413; 
ALEPH Collaboration, D. Buskulic et al, Z. Phys. C69 (1996) 393; 
DELPHI Collaboration, P. Abreu et al, Z. Phys. C68 (1995) 353; 
L3 Collaboration, M. Acciarri et al, Phys. Lett. B345 (1995) 589. 

[3] M. Anselmino et al., Quark Fragmentation into Vector and Pseudoscalar Mesons at LEP, 
CERN-PPE/97-136, submitted to Phys. Lett. B. 

[4] ALEPH Collaboration, D. Buskulic et al., Z. Phys. C62 (1994) 1; 

DELPHI Collaboration, P. Abreu et al, Z. Phys. C59 (1993) 533, erratum-ibid. C65 
(1995) 709. 

[5] F. Fiedler, Measurements of Hadronisation Aspects of Heavy Quark Production in Z° De- 
cays, to be published in the proceedings of the IVth International Workshop on Progress 
in Heavy Quark Physics, Rostock, Germany, September 20-22, 1997. 

[6] Yi-Jin Pei, Z. Phys. C72 (1996) 39. 

[7] L3 Collaboration, M. Acciarri et al., Phys. Lett. B396 (1997) 327. 
[8] OPAL Collaboration, G. Alexander et al., Z. Phys. C72 (1996) 1. 

[9] OPAL Collaboration, K. Ackerstaff et al., Measurement of /(c ^ D*+X), /(b -> D*+X) 
and r c c/r h ad using D*+ Mesons, CERN-PPE/97-093, accepted by Z. Phys. C. 

[10] OPAL Collaboration, K. Ackerstaff et al., Production of P-Wave Charm and Charm- 
Strange Mesons in Hadromc Z° Decays , CERN-PPE/97-035, accepted by Z. Phys. C. 

[11] OPAL Collaboration, P.P. Allport et al, Nucl. Instr. and Meth. A346 (1994) 476; 
OPAL Collaboration, P.P. Allport et al., Nucl. Instr. and Meth. A324 (1993) 34; 
O. Biebel et al., Nucl. Instr. and Meth. A323 (1992) 169; 
M. Hauschild et al., Nucl. Instr. and Meth. A314 (1992) 74; 
OPAL Collaboration, K. Ahmet et al, Nucl. Instr. and Meth. A305 (1991) 275. 

[12] OPAL Collaboration, G. Alexander et al., Z. Phys. C52 (1991) 175. 

[13] T. Sj6strand, Comp. Phys. Comm. 82 (1994) 74. 

[14] OPAL Collaboration, G. Alexander et al., Z. Phys. C69 (1996) 543. 

[15] C. Peterson, D. Schlatter, I. Schmitt and P. M. Zerwas, Phys. Rev. D27 (1983) 105. 

[16] J. Allison et al., Nucl. Instr. and Meth. A317 (1992) 47. 



25 



[17] M. Hauschild et al., Nucl. Instr. and Meth. A314 (1992) 74. 
[18] OPAL Collaboration, R. Akers et al, Z. Phys. C66 (1995), 19. 
[19] OPAL Collaboration, G. Alexander et al., Z. Phys. C73 (1997) 587. 
[20] Particle Data Group, Phys. Rev. D54 (1996) 1. 

[21] CLEO Collaboration, L. Gibbons et al., The Inclusive Decays B — > DX and B — > D*X 7 
CLNS 96/1454, hep-ex/9703006 (1996). 

[22] The LEP Collaborations ALEPH, DELPHI, L3, and OPAL, and the LEP Elcctroweak 
working group, Combining Heavy Flavour Electroweak Measurements at LEP, CERN- 
PPE/96-017 (1996). 

[23] OPAL Collaboration, K. Ackerstaff et al., Z. Phys. C74 (1997) 437. 

[24] The LEP Collaborations ALEPH, DELPHI, L3, and OPAL, the LEP Electroweak working 
group, and the SLD Heavy Flavour Group, A Combination of Preliminary Electroweak 
Measurements and Constraints on the Standard Model, CERN-PPE/96-183 (1996). 

[25] CLEO Collaboration, J. Gronberg et al., Phys. Rev. Lett. 75 (1995) 3232. 

[26] CLEO Collaboration, R. Balest et al., Phys. Rev. Lett. 72 (1994) 2328. 

[27] S. Godfrey and R. Kokoski, Phys. Rev. D43 (1991) 1679. 

[28] CLEO Collaboration, Y. Kubota et al., Phys. Rev. Lett. 72 (1994) 1972. 

[29] A. F. Falk and M. E. Peskin, Phys. Rev. D49 (1994) 3320. 

[30] F. Becattini, Universality of thermal hadron production in pp 7 pp and e + e~ collisions, hep- 
ph/9701275, v2, to be published in the proceedings of the XXXIII Eloisatron Workshop 
on Universality features in multihadron production and the leading effect, October 19-25 
1996, Erice, Italy. 



26 



A Model Independent Calculation of pP™ 

In this appendix, details of the calculations that lead to the dependence of the model indepen- 
dent value of Py im on the production of the broad L—l charmed mesons (equation (|T§D) are 
provided. 

As discussed in the text, the contributions from decays of D** mesons have to be taken into 
account when calculating P v nm . Two assumptions are made on D** production: 

• Since D** + mesons have not yet been observed in Z° — > cc decays, isospin invariance is 
assumed to be valid which yields equal production rates for corresponding D**° and D** + 
spin states. 

• For each light quark flavour q=u,d,s, two broad and two narrow L=l cq mesons are 
predicted. For the cu system, both narrow resonances (D? and Djf) have been measured 
at OPAL, whereas only one narrow cs resonance (D^) has been measured. Thus, the 
additional assumption is made that the relative production of the different L=l spin 
states does not depend on the flavour of the light quark. 

Under these two assumptions, the production of any L=l charmed meson can be expressed in 
terms of measured rates and two unknown parameters, /(c^Dq ) and /(c^D'j ). 

The OPAL production measurement |10| assumes a 100% branching for D£ -> D«K 



decays (i. e. to non-strange, charmed non-orbitally excited mesons). Therefore, any contribution 
from D** + — > DW + X decays is implicitly taken into account when using the measured 
multiplicity for an evaluation of P v nm . 

The quantity P v nm is calculated as follows: 

T/prim 

Py Tim = ~ (22) 

v ^/prim _|_ ^>prim v * 

+ P)eff _ ( T(s) ^U or P) ' 

where in the above formula, V and P stand for non-strange charmed vector and pseudoscalar 
mesons, respectively, whereas T( s ) denotes the sum of non-strange and strange charmed tensor 
mesons (D** and D**). Thus, 

im 2/(c^D*+) - 2/(c^D"°^D*°/+) - /(c^D b **+^D*°/+) 

v /(c->D°) +/(c^D+) -2/(c^D**°) -/(c^D**+) ' 1 ] 

where the factors of 2 are introduced from the assumptions f(c — > D*°) = f(c — ► D* + ) and 
/(c^D**°) = /(c^D**+). 

The symbol /(c^D** — >D* 0//+ ) is a shorthand notation for the hadronisation fraction of a 
charmed quark to a D** meson times the fraction of D** mesons decaying to a D*°/ + . Spin and 
parity conservation restrict the possible decays of the different D** spin states, as is illustrated 
in table || where a summary of predictions and measurements |20| for the four D**° spin states 



is given. Thus, the fraction of D** mesons decaying to a D* 0//+ depends both on the relative 
production of the different D** spin states and on the branching ratio BriJ)^ ~^ D*°/ + ), since 
for the D2 resonance, decays to both D* and D are allowed. Due to phase-space effects, the 
branching fraction Br** = Pr(D* 2 h — > D*°/ + ) of D*2~ mesons to a D*° or D* + is expected to 
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differ from the corresponding quantity Br** = Br(D*P — > D*°/ + ) for non-strange D** mesons, 
see table || 



name 


^0 




D° 




j = s q + L 


(1/2)+ 


(3/2)+ 


decay 


S-wave 


D-wave 


spin-parity J p 


0+ 


1+ 


1+ 


2+ 


decay channels 


Dtt 


D*7T 


D*7T 


Dtt, D*tt 


width (MeV) 


~ 100 


~ 100 


18.9±fj 


23 ±5 


mass (MeV) [0 


(not observed) 


2422.2 ±1.8 


2458.9 ± 2.0 



Table 6: Properties of the neutral, non-strange excited D mesons. For the decay channels, 
the same restrictions as shown here apply to D** + and D** + mesons, where in the case of D** + 
decays, pions have to be replaced by kaons. 



Using the shorthand notation /d x = f(c — > Dx), equation (|19|) can then be derived as 
follows: 



ppnm 



2/d. - 2 (/ Df + /„; + Br**f Df ) - ^ (/ Df + / D? + Br* s *f Df 

/d« + /d+ - 2 (/ D *o + / D /o + / D o + / D »o) - (/ D ,o + / D /o + / D o + / D *o 

-4 - -Rc/d'O 



(25) 



B - R c (/of + /d^ 

where in the transition to the second line, both numerator and denominator are multiplied by 
a factor of R c / (2 + j^E^boj' J • T ne parameters A and B are then given by 

2/(c^D* 0/+ ) - 2 (/( C ->D?) + /(c^D;°^D* 0/+ )) 

, /(c->p s " + : 

/(c^D**°) 

7(c-D°) + /(c-D+) - 2 f/(c-D°) + /(c^D* )' 



.4 



5 



respectively. 



R r 



/(c-d s +) + /(c^d: 2 + ^d* o / + ));/ 
/(c-d: 



and 



(26) 



/(c^D* 



/(c^D+)±/(c^D 



s2 



/ 



/(c^pr + ) 

/(c->D**°) 



(27) 
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Figure 1: The distributions of the four D*° selection variables m^o, cos#*, E 1 and E n o are 
shown for all preselected candidates. In each case, points with error bars correspond to the 
data and the open histogram to the simulation, scaled to the same number of entries. In figures 
(a) and (b), the contribution from D*° candidates with a correctly reconstructed D°-^K~7r + 
decay is shown as the hatched area. The cross-hatched areas in figures (c) and (d) correspond 
to 10 times the contribution of correctly reconstructed photons or ir° mesons from a D*° decay. 
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Figure 2: The distributions of the likelihood functions C for the two D*° decay channels. 
Points with error bars represent the data and open histograms the simulation, scaled to the 
same number of entries. The hatched histograms show candidates in the simulation with a 
correctly identified photon or 7r° from a D*° decay; for candidates entering the cross-hatched 
histogram, the D° — > K~n + decay is also correctly reconstructed. The vertical lines represent 
the cuts on the likelihood functions. No scale is given for the likelihood values since only the 
relative position of the cut is relevant. 
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Figure 3: The Am(D°,7) and Am(D°,7r°) mass difference distributions in the data and the 
simulation. 

In (a), (b), (e), and (f), the distributions obtained in the signal selection procedure are shown as 
solid histograms, and those obtained in the background selection procedure as points with error 
bars, where the relative normalisation of the latter distribution has been determined in the fit. 
In addition, the hatched histograms show the distributions of signal candidates reconstructed 
in the simulation with the signal selection procedure. 

In (c), (d), (g), and (h), the corresponding background subtracted distributions are shown 
together with the fit results; the error bars show the statistical errors only. 
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Figure 4: The mass difference distribution for the decay D* + — > D+7. As before, in (a) and 
(b), the solid histogram corresponds to the signal sample, and points with error bars to the 
background sample. The signal contribution in the simulation is indicated by the hatched his- 
togram. 

The corresponding background subtracted distributions are shown in (c) and (d) together with 
the fitted signal parametrisation. The shaded Gaussian in figure (c) shows the expected con- 
tribution from D* + — >D+7r° decays which is fixed in the fit. 
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Figure 5: The dependence of Py im on the two unmeasured multiplicities R c x /(c- 
R c x /(c— >D' 1 °) is shown in figure (a). Each of these quantities is varied between and twice 
the value as expected from a spin counting picture. Plots (b) and (c) show the dependence on 
each one of these multiplicities when the other is fixed at the spin counting expectation. The 
error contours correspond to one standard deviation. They include the experimental error on 
the relative production of D 1 ^ and D^ mesons as measured in flIQfl . 
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